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Abstract. The heat capacity anomaly at the transition to superconductivity of the layered 
superconductor MgB 2 is compared to first-principles calculations with the Coulomb repulsion, 
/J,*, as the only parameter which is fixed to give the measured Tc. We solve the Eliashberg 
equations for both an isotropic one-band and a two-band model with different superconducting 
gaps on the vr-band and cr-band Fermi surfaces. The agreement with experiments is 
considerably better for the two-band model than for the one-band model. 
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1. Introduction 

The nature of the superconducting state in MgB2 has been characterized by a broad range of 
experimental and theoretical methods and many basic properties have been unambiguously 
established since the discovery of the 40 K superconductor MgB2 by Nagamatsu and 
collaborators 

While electron-phonon coupling as the underlying pairing mechanism has been 
pinpointed by a large B isotope effect on Tc proving B related vibrations to be essential 
[0, 1^] the nature of the order-parameter (viz. the superconducting gap) has remained a 
matter of debate. The order-parameter has been intensively investigated by tunnelling and 
point contact spectroscopy [P]-p^ as well as by high-resolution photoelectron spectroscopy 
[|T7|,|rS]]. While these techniques show an energy gap in the quasiparticle spectrum most 
likely of s-wave symmetry the magnitude of the gap, A(0), itself remained an open 
question: tunnelling experiments initially revealed a distribution of energy gaps with lower 
boundary 2Ai(0)/fcBT'c ~ 1.1 and upper boundary 2A2(0)//cb^c ~4.5. These values 
are either considerably lower or distinctly larger than the weak coupling BCS value of 
2A(0)//cbTc = 3.53 and these controversial findings have been discussed in terms of gap 
anisotropy or more recently attributed to the presence of two gaps or multiple gaps [[T2|,[r8|]. 
The analysis of the electronic Raman continuum of MgB2 by Chen et al. [|T3]] also pointed 
to the presence of two gaps with gap values within the limits indicated by the tunnelling 



experiments [[T9|]. 

While these experiments employ surface sensitive techniques to determine the gap 
properties, evidence for multigap behaviour emerges also from methods like heat capacity or 



fxSR measurements probing true bulk properties [ pO| ] . In the early heat capacity experiments 
the typical jump-like anomaly is seen at the magnitude ACp/T^ of which amounts at best 
to only about 70-80% of the value 1.43 7^ (Tc) predicted by weak-coupling BCS theory 
P,p]-P5|]. 7^(Tc) is the Sommerfeld constant in the normal state which was obtained from 
heat capacity measurements in high magnetic fields and which was determined to be 2.7- 
3 mJ mol^^ K^^. The shape of the heat capacity anomaly compares reasonably well with 
BCS-type behaviour assuming 2A(0)/A;b^c = 3.53 with appropriately adjusted magnitude. 
An improved fit of the detailed temperature dependence of the heat capacity anomaly was 
obtained when calculating the heat capacity within the a-model assuming a BCS 
temperature dependence of the gap but with an increased ratio 2A(0)//cb^c = 4.2(2) [^. 
This result matches very well with the upper limit of the gap value consistently found in 
the tunnelling experiments and was suggested as an evidence that MgB2 is in the moderately 
strong coupling limit. More recently, the excess heat capacity observed close to TJA by 
Bouquet et al. [ p2[ ] and Wang et al . [ ^4] ] has been attributed to a second smaller gap. Fits 
with a phenomenological two-gap model assuming that the heat capacity of MgB2 can be 
composed as a sum of the two individual heat capacities gave very good description with 



gap values of 2Ai(0)//cb^c = 1-2(1) and 2A2(0)/A;bTc ^ 4 p7|]. Recent muon-spin-relaxation 



measurements of the magnetic penetration depth are consistent with a two-gap model [ P8| ] 



Theoretically multigap superconductivity in MgB2 was first proposed by Shulga et al. 
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to explain the behaviour of the upper critical magnetic field [gyp. Based on the electronic 
structure the existence of multiple gaps has been suggested by Liu et al. in order to explain 
the magnitude of Tc [PPj]. The electronic structure of MgB2 contains four Fermi surface sheets 
[PT|]. Two of them with 2D character emerging from bonding a bands form small cylindrical 
Fermi surfaces around F - A. The other two originating from bonding and antibonding tt bands 
have 3D character and form a tubular network. Liu et al. from first-principles calculations 
of the electron-phonon coupling conclude that the superconducting gap is different for the 
individual sheets and they obtain two different order parameters, a larger one on the 2D Fermi 
surface sheets and a second gap on the 3D Fermi surfaces, the latter was estimated to be 



approximately a factor of three reduced compared to the former [ pOQ . 

In the present paper we calculate the specific heat capacity from the spectral Eliashberg 
function a^{ijj)F{uj) first in the one-band model using the isotropic a'^{ijj)F{ijj) as given by 
Kong et al. [ |32[ ] . Then we calculate the heat capacity in a two-band model by reducing the 
16 Eliashberg functions al-{uj)Fij{uj) appropriate for the four Fermi surface sheets into four 
Eliashberg functions corresponding to an effective-two-band model with a cr-band and vr-band 
only. From the solution of the Eliashberg equations we obtain a superconducting gap ratio 
Act/ — 2.63 in good agreement with the experimental data. The two-band model explains 
the reduced magnitude of the heat capacity anomaly at very well and also reproduces the 
experimental observed excess heat capacity at low temperatures. 

2. Theory 

2.1. One-Band Model 

First we discuss the specific heat in the isotropic single band model with a strong 
(intermediate) electron-phonon interaction (EPI). In the normal state and in the adiabatic 
approximation the electronic contribution to the specific heat is determined from the 
Eliashberg function a'^{ijj)F{ijj) by the expression 

Cfj{T) = (2/3)7r2iV(0)4r (1) 

j'OO 

l + {<6/TxkBT) I f{uj/2TikBT)a'^{uj)F{uj)duj 



where A^(0) is a bare density of states per spin at the Fermi energy. The kernel f{x) is 
expressed in terms of the derivatives of the digamma function ^(x) 

f{x) = —x — 2x'^lm^'{ix) — x"^ Re ^" (ix) . (2) 

At low temperatures the specific heat has the well known asymptotic form: C^{T 
0) = (1 + A)7oT, where A = 2 du!U!~^a'^{u!)F{uj) is the electron-phonon coupling 
constant, and 70 = 27r^A;|A^(0)/3 is the specific heat coefficient for noninteracting electrons. 
At higher temperatures the specific heat differs from this trivial expression (see, the discussion 
in reference [p4|]). 

In the superconducting state an expression for the specific heat obtained by Bardeen 
and Stephen [551] which is based on an approximate sum rule has often been used. We shall 
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however use an exact expression for the thermodynamical potential in the electron-phonon 
system which is based on the integration of the electronic Green's function over the coupling 
constant 

1 



t{x)G{x) 



(3) 



where x is dimensionless, G{x) and S(a;) are the exact electron Green's function and the 
self-energy, respectively, for a coupling constant of x ■ A. The functions Vt^^^ and Vt^^^ 
are the thermodynamic potentials for noninteracting electrons and noninteracting phonons, 
respectively. Some further arithmetics leads to the expression for the difference in free 
energies, Fn and Fs, of the superconducting and normal state [B^ 



ttN{0)T 



E 



7N, 



11 



2a;2[(zS(c^„))'-l]+y2 



(4) 



where Z{ujn) is a normalization factor, ipn = A„/Z(co'„) is an order parameter, and A„ is the 
gap function. 

The specific heat at temperature, T, is calculated according to: 

ACei(T) = Td\FM - Fs)/dT\ (5) 

The specific heat jump ACci(Tc) at T = Tc is determined by the coefficient (3 = 

TAC,i{T,)/2 in Fn-Fs = (3t^ where t = (T, - T)/T,. 



2.2. Two-Band Model 

We have calculated the 16 Eliashberg functions afj{uj)Fij(uj) where i and j label the four 
Fermi surface sheets and thereafter combined them into four corresponding to an effective 
two-band model which contains only a a- and vr-band. Their respective densities of states 
at the Fermi energy have values of A^<^(0) = 0.300 states/cell-eV and iV^(O) = 0.410 
states/cell-eV. Similar coupling constants X^^^X^^^^X^^, and Xt,„ which are required for a 



two-band model were calculated earlier in reference [pO|]. The procedure of reducing the 
16 Eliashberg functions of the real 4 band system due to the 4 different Fermi surface sheets 
to an effective two-band model with only four coupling constants Xij is an approximation 
which is based on the similarity of the two cylindrical and the two three-dimensional sheets 
of the Fermi surface requiring the same physical properties in both cr-bands or both tt -bands. 
More details can be found elsewhere p7| ] . 

The four Eliashberg functions alj{uj)Fij{ijj) for the effective two-band model are shown 
in figure [1]. The most significant contribution comes from the coupling of the bond stretching 
phonon modes to the cr-band. The coupling constants corresponding to the superconducting 
Eliashberg functions have been calculated to be: X^a = 1.017, Ajr^ = 0.448, A(j7r = 0.213, 
and A^o^ = 0.155. The small difference to the values given in reference [ PD| ] may be attributed 
to the different first-principles methods used in the calculation of the Eliashberg functions. 
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Figure 1. The four superconducting Eliashberg functions F{uj) obtained from first- 
principles calculations for the effective two-band model and the isotropic Eliashberg function 
for the one-band model. The coupling constant of the isotropic one-band model has a value of 



Besides the spectral functions we need to know the the Coulomb matrix element 
Hij. With the help of the wavefunctions from our first-principles calculations we can 
approximately calculate the ratios for the yU-matrix [|J7[]. The aa-, nn- and avr-values were 
in the ratio 2.23/2.48/1. This allows one to express f^*j{ujc) by these ratios and one single 
free parameter which is fixed to get the experimental Tc of 39.4 K from the solution of 
the Eliashberg equations. The ij*{uJc) matrix elements determined by this procedure are 
/i:,(cue)=0.210, /i:^(c^,)=0.095, /i;,(cue)=0.069, and K^(c^,)=0.172. 

Using our calculated Eliashberg functions on the imaginary (Matsubara) axis together 
with the above matrix iJ*ij{uJc) we obtain the gap values A^r = lini Ao^(i7rT) ~ 7.1 meV, and 
~ 2.7meV, which corresponds to 2Ao-/A;bTc=4.18 and 2A^/A;bTc=1.59. The temperature 
dependence of the superconducting gaps is shown in figure 0. The filled circles (squares) 
display the gap for the 2D a- (3D tt-) band. Due to interband coupling between the bands 
both gaps close at the same critical temperature. For a comparison also the BCS curve (line) 
is shown for a single gap (one-band model) which closes at Tc=39.4 K. The corresponding 
single BCS gap would be 6 meV. 

The extension of equation |T] to the two-band model gives 

Cl{T) = '^NtotmlT + 47rA;B[iV.(0)(/^. + h^) + iV^(0)(/^^ + /^.)](6) 

where = f{iu/2nkBT)afj{uj)Fij{iu)duj = tt, a), and the function f{x) is given by 
equation |^. 

The generalization of the superconducting free energy (^ to the two band model is 
straightforward and the heat capacity was obtained according to equation |5|. 



Specific heat of MgB2 from first-principles calculations 



6 




Figure 2. The temperature dependence of the the superconducting gaps from the solution 
of the two-band EHashberg equations. The values of the gaps at T=0 K were obtained 
as A,, (T = 0) = 7.1 meV and A„ (T = 0) = 2.7 meV. The BCS value for the gap that 
con-esponds to Tc=39A K is 6.0 meV at K. 



3. Comparison with Experiment 

For the comparison with experiment we have selected data obtained by our group and 



by Bouquet et al . p2|]. The anomaly clearly visible at in the zero-field data is suppressed 
by a magnetic field of 9 Tesla in both experiments. In figure ^ we display the difference 
ACp = Cj,(OTesla) — Cp(9Tesla). The anomalies at detected by both groups clearly have 



a different magnitude, the one described in reference p2|] amounts to 133 mJ/mol K at T, 



and represents the largest specific heat capacity anomaly reported for MgB2 so far [p8|]. The 
ACp{Tc) reported by our group is somewhat smaller, however, the shape of the anomalies 
close to Tc is very similar for both samples. In fact, fitting the anomalies with the «-model 
revealed an identical ratio 2A//cbTc=4.2 with A=7 meV for both samples pT|,p2|,^. 



First we will try to discuss the experimental results in terms of a conventional one-band 
model. The specific heat in MgB2 was calculated using the isotropic spectral Eliashberg 
function a'^(uj)F(Lu) of Kong et al. [B^. This function yields an electron-phonon coupling 



constant A = 0.87 and together with a Coulomb pseudopotential of jj,* = 0.1 yields Tc = 40 
K. The calculated specific heat at is -f^{Tc) = 1.947o = 3.24 mJ/molK^ with 7o=1.67 
mJ/molK^ from the band structure calculations of reference [|3l|,^^. The specific heat jump 
at Tc equals AC ~ 196 mJ/mol K, which is a factor of 1.5-2 larger than the experimental 
values It corresponds to AC/(7^ (TjTc) ~ 1.51 compared to the BCS value of 

1.43. The difference ACei{T) = C|(r) - C^Y(T) is shown in figure | (dashed-dot line) 
in comparison with the experimental data. Not only the size of the jump disagrees with the 
experiment, but also the behaviour at low temperatures is different. The latter is connected 
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Figure 3. Experimental data of the heat capacity difference ACp — Cp (OTesla) — Cp (9Tesla) 
from reference [ pl| ] (o) and from reference (A). The dashed line is the theoretical result 
of the one-band model and the thick solid line corresponds to the two-band model from the 
solution of the Eliashberg equations. The two-band model reproduces much better the specific 
heat jump as well as the low temperature behaviour. 



with the fact that at low temperature equation (Q) for a single band model leads to the standard 
exponential dependence ~ T^^^"^ exp(— A/T), while the experimental data show a more 
complicated behaviour. Clearly there exists a discrepancy between experimental data and a 
theoretical one-band model. 

The solid line in figure ^ represents the theoretical results for the two-band model as 
described above. The low temperature behaviour is in much better agreement with the 
experiment. The specific heat jump is now significantly reduced in comparison with a single 
band model and reproduces surprisingly well the experimental data of reference [|2^]. With 
the data given above we obtain from our theoretical calculation an electronic heat capacity in 
the normal state of 7^(T'c) = C^i(Tc)/Tc ~ 3.24 mJ/mol K^, the same value as for the one- 
band model. The absolute value of the specific heat jump in the two-band model is AC ~ 
125 mJ/mol K, corresponding to AC/(7^ (Tc)Tc) ^ 0.98 which is now smaller than the BCS 
value. 

We would like to emphasize here that no fitting is involved in the theoretical calculations. 
The only free parameter which is in the Coulomb matrix elements is already determined by 
the experimental Tc of 39.4 K. 

One could expect that the difference between the theoretical results of the effective two- 
band model and our experimental data may be attributable to a different amount of impurities 
in the samples compared to the samples of reference In the one-band model the critical 
temperature Tc as well as the value and the temperature dependence of ACp(T) are not 
affected by non-magnetic impurities (Anderson theorem). This is in complete contrast to the 
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situation for the two-band model, where both quantities are strongly dependent on interband 
impurity scattering. Interband impurity scattering leads to averaging of the gaps and thus to 
the increase of ACp/'j^ (Tc)Tc ratio. On the other hand, due to decrease of Tc, the specific heat 
jump only depends weakly on the scattering strength. In order to investigate the dependence 
of Tc and of ACp on the interband impurity scattering we included the effect of interband 
impurities in the Eliashberg equations. The results show that even for rather strong impurity 
scattering l/2r = SttTcq — 371 K, which leads to a drastic change of the critical temperature 
(decreasing to Tc=29.4 K) and strong averaging of the gaps, the specific heat jump remains 
practically unchanged ACp ^ 120 mJ/mol K. This corresponds to a ratio ACp/'y{Tc)Tc ~ 
1.48, which is very close to the corresponding value of a single gap model. Therefore, 
interband impurity scattering can explain the change of Tc in different samples, but is not 
responsible for the observed different values of the specific heat capacity anomaly at Tc. 

We have shown that a complete theoretical calculation from first-principles using an 
effective two-band model can explain the major features in the specific heat measurement 
of MgB2 surprisingly well. The presented theoretical framework goes beyond a simple 
phenomenological two-gap model because interband effects are included explicitely and 
no fitting to experimental results has been performed. The reduced value of the heat 
capacity anomaly at Tc as well as the low temperature behaviour are in excellent agreement 
with experimental results. The same first-principles approach using exactly the same 
Eliashberg functions and Coulomb matrix elements has been used in order to explain optical 
measurements [ |J7| ] and tunnelling experiments of the interesting superconductor MgB2. 
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